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ABSTRACT 


'-y 

An  analytical/computational  technique  has  been  developed  for  determining 
the  geometrically  and  constitutively  nonlinear  response  of  an  infinite 
cylindrical  shell  to  a  transverse,  transient  acoustic  wave.  Shell  be¬ 
havior  has  been  treated  through  utilization  of  the  nonlinear  structural 
analyzer  DYNAPLAS  II,  while  the  fluid-structure  interaction  has  been 
treated  in  accordance  with  both  the  exact  residual  potential  formulation 
and  the  doubly  asymptotic  approximation.  Numerical  results  produced 
through  application  of  the  approximation  differ  significantly  from  the 
corresponding  exact  results. 
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Section  1 


INTRODUCTION 

About  ten  years  ago,  complete  analytical  solutions  for  the  two-dimen¬ 
sional,  plane-strain  response  of  a  linear-elastic,  circular  cylindrical  shell  to 
a  transient  acoustic  wave  first  appeared  [1,2] .  The  solution  techniques  Involved 
Fourier  decomposition  with  respect  to  the  circumferential  coordinate,  followed 
by  the  Introduction  of  differential  and  integral  equations  to  treat  the  fluid- 
structure  interaction  on  a  harmonic-by-harmonic  basis.  Recently,  these  solu¬ 
tions  were  extended  to  include  geometrically  nonlinear  behavior  of  the  shell 
in  order  to  examine  the  elastic  dynamic  instability  of  submerged  cylindrical 
shells  [3],  The  purpose  of  the  present  study  has  been  to  extend  the  solutions 
still  further  in  order  to  examine  inelastic  shell  response. 

Analytical  solutions  to  idealised  problems  offer  two  primary  benefits  to 
the  technical  community.  First,  they  facilitate  the  development  of  physical 
insight  into  the  phenomenology  involved.  This  is  not  only  because  of  the 
relative  simplicity  of  the  problems  treated,  but  also  because  of  the  opportun¬ 
ity  to  perform  extensive  parameter  studies  at  modest  cost.  Second,  analytical 
solutions  serve  to  provide  check  problems  for  computer  codes  developed  to  treat 
much  more  complex  problems.  Clearly,  no  code  should  be  applied  to  complex 
problems  before  it  is  thoroughly  tested  on  a  number  of  check  problems. 

The  present  study  might  not  be  considered  by  some  to  be  an  analytical 
study,  in  that  the  structural  analysis  code  DYNAPLAS  [4,5]  is  used  to  treat  the 
geometrically  and  constitutively  nonlinear  behavior  of  the  shell.  DYNAPLAS 
was  selected  for  this  purpose  because  it  is  well  established  and  it  is  based 
upon  circumferential  Fourier  decomposition  of  shell  response. 

As  discussed  in  Section  2,  "live-load"  terms  have  been  appended  to  the 
DYNAPLAS  equations  of  motion  to  account  for  geometrically  nonlinear  loading 
effects.  These  terms,  which  are  derived  in  [3],  are  treated  as  pseudo-forces 
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during  Che  solution  process.  Also  discussed  in  Section  2  is  the  approach  se¬ 
lected  for  treatment  of  the  fluid-structure  interaction,  viz.,  the  "residual 
potential  formulation"  (RPF)  of  [1].  This  approach,  which  is  based  upon  cir¬ 
cumferential  Fourier  decomposition,  constitutes  an  exact  treatment  of  fluid- 
structure  interaction  for  a  surrounding  acoustic  medium. 

Numerical  results  are  presented  in  Section  3  for  two  check  problems  de¬ 
signed  to  verify  the  solution  procedure  and  an  inelastic  response  problem 
designed  to  highlight  the  phenomenology  of  Interest.  The  response  problem 
involves  a  shell  with  elastic/perfectly  plastic  material  behavior  that  is 
characterized  by  a  static,  elastic,  critical-buckling  pressure  and  an  axi- 
symmetric  elastic-limit  pressure  that  are  virtually  equal.  The  shell  is  ex¬ 
cited  by  a  rectangular  acoustic  wave  with  a  pressure-magnitude  equal  to  four 
times  the  shell's  elastic  critical-buckling  pressure  and  a  width  equal  to 
four  shell  radii. 

Numerical  results  are  also  presented  in  Section  3  for  the  inelastic  re¬ 
sponse  problem  with  the  fluid-structure  interaction  treated  in  accordance  with 
the  "doubly  asymptotic  approximation"  (DAA)  [6,7].  This  approximation,  which 
is  the  basis  for  fluid-structure  interaction  analysis  in  a  number  of  existing 
codes  [8-13],  is  asymptotically  exact  for  both  low-  and  high-frequency  fluid 
motions,  effecting  a  smooth  transition  in  the  intermediate  frequency  range. 

Its  computational  advantage  is  that  it  may  be  expressed  as  a  matrix  ordinary 
differential  equation  without  requiring  discretization  of  the  infinite  volume 
of  fluid  surrounding  the  structure. 

Section  4  completes  the  report  with  a  summary  of  the  work  performed  and 
a  list  of  conclusions. 
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Section  2 

GOVERNING  EQUATIONS 

Consider  the  two-dimensional,  plane-strain  motions  of  the  submerged, 
infinite,  circular  cylindrical  shell  of  Figure  1.  The  shell  is  excited  by  a 
transient  acoustic  wave  that  first  contacts  the  shell  at  8  *  tt.  During  the 
resulting  fluid-structure  interaction,  shell  behavior  may  involve  both  geo¬ 
metric  and  constitutive  nonlinearity. 

2.1  STRUCTURAL  EQUATIONS 


The  decomposition  of  shell  and  fluid  response  into  circumferential 
Fourier  harmonics  yields 
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where  $(r,9,t)  is  the  velocity  potential  for  the  acoustic  field.  Under  such  a 
decomposition,  the  structural  equations  of  motion  for  the  axisymmetric  and  nth 
nonaxisymmetric  response  harmonics  are  [3,4] 
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where  p  (r,t)  is  the  nth  harmonic  of  fluid  pressure  and  f  ,  f _ and  f_  are 
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stiffness-force  harmonics  computed  within  DYNAPLAS  that  involve  linear-elastic 
geometrically  nonlinear  and  constitutively  nonlinear  behavior. 

The  terms  on  the  right  sides  of  (2)  that  involve  pressure-displacement 
products  are  "live-load"  terms  appropriate  to  moderate  displacements  and 
rotations  of  the  shell  surface.  They  account  for  the  major  effect  of  finite 
membrane  strain  and  finite  rotation  of  shell  normals  on  shell  excitation  by 
a  normal  pressure  loading,  as  well  as  the  primary  effect  of  finite  shell  dis¬ 
placement  on  the  loading  produced  by  the  spatially  varying  acoustic  field 
[3].  Without  these  terms,  the  right  sides  of  (2)  would  involve  only  "dead 
load"  terms  appropriate  to  infinitesimal  shell  motions.  As  DYNAPLAS  does  not 
possess  live- load  capability,  the  live-load  terms  have  been  treated  as  pseudo¬ 
forces  generated  by  auxiliary  software. 


2.2  FLUID  STRUCTURE  INTERACTION 


A  residual  potential  formulation  (3? F)  of  the  fluid-structure  interac¬ 
tion,  which  constitutes  an  exact  treatment,  may  be  constructed  as  follows  [1]. 
First,  fluid  pressure  and  radial  fluid-particle  velocity  are  expressed  as 
derivatives  of  the  fluid  velocity  potential  as 

p  *  -p4> 

(3) 

u  -  3$/3r 

Second,  the  total  acoustic  field  is  treated  as  the  superposition  of  the  acous¬ 
tic  field  for  the  known  incident  wave  and  the  acoustic  field  for  the  unknown 
scattered  wave,  i.e.. 


$(r,e,t)  -  bjCr.e.t)  +  4>s(r,0,t) 


(4) 


Third,  compatibility  of  radial  fluid-particle  velocity  and  radial  shell  vel 
ocity  is  enforced  at  the  wet  surface  of  the  shell  as 


u(a,e,t) 
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Finally,  Che  wave  equation  and  radiation  condition  for  each  circumferential 
harmonic  of  the  scattered  wave  are  replaced  by  the  equivalent  residual-poten¬ 
tial  equation 
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in  which  the  residual  potential  is  given  by  the  convolution  relation 
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where  the  r  are  characteristic  functions  that  resemble  step-exponential 
tl 

functions  tl]. 


Equations  (3)  -  (7)  may  be  utilized  to  produce,  for  each  circumferential 
harmonic,  the  fluid-structure- interaction  relations 
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where  *  A_  (a,t)  is  obtained  from  (7). 

Kn  Kn 


2.3  SOLUTION  PROCEDURE 


In  principle,  (2),  (7)  and  (8)  may  be  solved  simultaneously  in  direct 
fashion  for  a  selected  number  of  Fourier  harmonics;  Fourier  superposition  in 
accordance  with  (1)  then  produces  response  histories  of  interest.  In  prac- 
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tice,  however,  care  must  be  taken  to  provide  an  implmentation  that  yields 
solutions  of  satisfactory  accuracy  at  acceptable  computational  cost.  As 
verified  in  Section  3,  such  an  implemnetation  has  been  constructed  in  non- 
dimensional  form  in  accordance  with  the  normalization 


^  a  2  » 

p  *  p  /pc  ,  w  •  w/a,  t  *  ct/a 


The  nondimensional  response  equations  underlying  the  implementation 
used  for  computations  are  [cf.  (2),  (7)  and  (8)] 

(— )  (~)w  +  w  +  f  (v,w)  =*  +  u_  -  -i-  +4_ 

v  p  '  va'  o  o  o  *  *Io  Io  2  *So  *Ro 


(“)  +  w  +  f  (v,w) 

p  &  n  n  wn 


►in  +  uIn  “  2  *Sn  +  *Rn 


+  ($T(>  +  )(nv  +  w  )  +  w  w  n  >  1 
io  So  n  n  on 


o  h  .  (10) 

(T)(a)vn  +  fvn(v*w)  “  (*Io  +  *So)(vn  +  «wn)  ’  n  5  1 


►sn  +  2  '►sn 


u_  -  &  +  <b. 

^n  n  I 


?Rn  ’  ”rn  *  ^Sn 


where  all  circumflexes  have  been  dropped,  wq  =  dwQ/dt,  etc.,  and  the  asterisk 

denotes  temporal  convolution,  as  defined  in  (7).  For  n  »  0,  (10)  contains 

three  equations  for  the  three  unknowns  w  ,  and  ,  while  for  each  n  5  1, 
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(10)  contains  four  equations  for  the  four  unknowns  w  ,  v  ,  and  .  From 
a  computational  standpoint,  the  primary  advantage  of  (10)  is  that  the  terms 
involving  unknowns  on  the  right  sides  of  the  first  three  of  (10)  are  more 
slowly  varying  than  the  acceleration  and  velocity  terms  on  the  left  sides. 


The  combined  response  equations  (10)  were  solved  by  step-by-step  numeri¬ 
cal  integration  as  follows.  The  first  three  of  (10)  were  solved  with  the  half- 
step  central-difference  algorithm  [14],  which  was  introduced  into  DYNAPLAS  for 
reasons  given  in  Section  3.  This  modification  also  involved  the  introduction 
of  structural  damping  analysis  capability  into  DYNAPLAS  to  treat  the  velocity 
terms  on  the  left  sides  of  the  first  two  of  (10).  The  fourth  of  (10)  was 
solved  with  a  fourth-order  Runge-Kutta  scheme  [15],  as  in  [3].  Finally,  the 
last  of  (10)  was  solved  with  trapezoidal  integration,  as  in  [3].  This  pro¬ 
cedure  possessed  satisfactory  accuracy,  stability  and  efficiency  characteris¬ 
tics  . 


Section  3 

NUMERICAL  RESULTS 


As  seen  in  Section  2,  the  derivation  of  the  structural  and  fluid-structure- 
interaction  equations  in  accordance  with  [3]  and  [4]  is  relatively  straightfor¬ 
ward.  The  establishment  of  a  satisfactory  solution  procedure,  which  involves 
the  interfacing  of  two  distinct  computer  codes,  is  not  straightforward.  Pre¬ 
sented  in  this  section  are  numerical  results  for  two  check  problems  designed  to 
verify  the  solution  procedure  and  a  response  problem  representative  of  the  phe¬ 
nomenology  of  interest,  viz.,  the  geometrically  and  constitutively  nonlinear 
behavior  of  shock-wave-excited,  submerged  shells.  All  results  are  presented 
in  nondimens ional  form  normalized  in  accordance  with  (9). 


3.1  CONTRIVED  DRY  CHECK  PROBLEM 


Consider  an  infinite  elastic,  circular  cylindrical  shell  responding  in 
plane  strain  to  the  prescribed  surface  pressure  loading 


p(0,t)  -  -  PQ  H(t)  H(to-t)  -  PQ  H(t-t0) 

o 

t_tl 

-P  — —  H(t-t. )  cos  n0,  n  >  2 
n  t  l 
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where  H(t)  is  the  Heaviside  step-function.  For  t^  >  tQ,  (11)  constitutes  an 
axisymmetric  ramp-step  loading  followed  by  a  nonaxisymmetric  ramp  loading  of 
specified  harmonic  index.  If  t^  >>  tQ  and  an  axisymmetric  damping  mechanism  is 
present  (which  is  the  situation  considered  here) ,  axisymmetric  shell  response 
reaches  its  static  asymptote  before  the  nonaxisymmetric  loading  is  applied. 
Subsequent  nonaxisymmetric  shell  response  then  constitutes  that  of  a  hydro¬ 
statically  pressurized  shell. 


As  discussed  in  [3],  nonaxisymmetric  shell  response  to  (11)  is  essentially 
inextensional.  Hence,  for  t^  >>  tQ  and  appreciable  damping  of  axisymmetric  shell 
response,  nonaxisymmetric  response  to  (11)  is  governed  by  the  equation  [3] 
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where  c  *  [E  /p  (1-v  )]^  is  the  plate  velocity  of  the  shell  material  and 
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P  »  (n  —  1) (p  /p)(c  /c)  (h  /12a  )  is  the  static  critical  buckling  pressure 
cn  o  o 

for  the  nth  circumferential  harmonic.  The  solution  to  (12)  for  quiescent 
initial  conditions  is 
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where  the  response  slope  Sn  and  the  natural  period  are  given  by 
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DYNAPLAS  II  and  DYNAPLAS  III  computations  were  performed  for  n  *  2-5  on 
the  basis  of  (2)  and  (11)  with  3pQ/3r  =  0;  the  parameter  values  chosen  for  the 
computations  are  shown  in  Table  1.  The  Houbolt  time- integration  scheme  was 
used,  which  is  characterized  by  high  algorithmic  damping  for  time-increment/ 
natural-period  ratios  greater  than  0.1  [16].  Hence,  with  the  time  increments 
of  Table  2  constituting  one-half  or  one-quarter  of  the  period  of  the  n  *  0 
shell  mode,  hydrostatic  conditions  were  reached  well  before  application  of  the 
nonaxisymmetric  loading  in  every  calculation. 

The  response-slope  and  natural-period  values  gleaned  from  the  response 
computations  are  shown  in  Table  2.  Good  agreement  is  observed  between  the 
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DYNAPLAS  II  and  analytical  values.  The  discrepancies  that  exist  are  primarily 
due  to  the  relatively  poor  accuracy  characteristics  of  the  Houbolt  scheme  [16]. 
Appreciable  discrepancies  exist  between- the  DYNAPLAS  III  and  analytical  values 
for  n  *  2,  however.  Furthermore,  no  DYNAPLAS  III  values  are  shown  for  n  5  3, 
because,  for  these  harmonics,  the  code  erroneously  produced  nonaxisymmetric  re¬ 
sponse  results  for  t  <  t^,  when  only  the  axisymmetric  loading  was  present. 

It  should  be  mentioned  that  a  "dead-load"  version  of  Table  2  was  also  com¬ 
piled,  corresponding  to  neglect  to  the  pressure-displacement  products  on  the 
right  sides  of  (2) .  This  was  done  to  ensure  that  no  discrepancy  between  a 
DYNAPLAS  result  and  an  analytical  result  could  be  attributed  to  the  auxiliary 
live- load  software  mentioned  in  Subsection  2.1.  The  dead-load  discrepancies 
observed  were  entirely  consistent  with  those  of  Table  2. 

In  view  of  Table  2,  DYNAPLAS  II  was  selected  as  the  structural  analyzer 
for  the  present  study.  This  was  originally  somewhat  of  a  disappointment, 
because  DYNAPLAS  II  does  not  treat  multilayer  shells,  whereas  DYNAPLAS  III 
does.  Multilayer  analysis  capability  was  desired  in  order  to  use  a  sandwich 
shell  as  a  plane-strain  model  for  a  stiffened  shell,  as  discussed  in  [1].  This 
modeling  problem  was  satisfactorily  overcome,  however,  as  described  below. 

Because  of  the  rather  severe  accuracy  limitations  of  the  Houbolt  scheme, 
the  half-step  central-difference  scheme  [14]  was  introduced  into  DYNAPLAS  II 
to  integrate  the  augmented  structural  response  equations  given  by  the  first 
three  of  (10).  This  produced  much  better  accuracy  relative  to  elastic-shell 
results  presented  in  [3]  than  that  provided  by  the  whole-step  central-difference 
scheme  already  existing  in  DYNAPLAS  II.  Hence  accuracy  requirements,  which 
mandated  augmentation  of  the  structural  response  equations,  also  mandated  re¬ 
structuring  of  the  DYNAPLAS  II  integrator. 

Finally,  one  additional  note  regarding  DYNAPLAS  calculation  should  be  men¬ 
tioned.  It  was  found  that  the  default  operation  involving  static  condensation 
of  the  rigid-body  shape- function  coefficients  (the  S's  of  p  24  in  [4])  produced 
an  erroneous  linear  stiffness  matrix  for  n*l.  The  entry  of  "1"  in  Column  30 
of  Card  i  in  the  SAMS0R4  input  deck  overrides  the  default  operation,  thereby 
correcting  the  problem. 
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3.2  REPRESENTATIVE  SHOCK-WAVE-EXCITATION  PROBLEM 

Consider  a  submerged,  infinite,  circular  cylindrical  shell  responding  in 
plane  strain  to  an  incident  wave  of  rectangular  pressure-profile.  The  desired 
shell  would  be  a  sandwich  shell,  exhibiting  the  enhanced  flexural  stiffness 
properties  characteristic  of  stiffened  shells  [1].  Because  DYNAPLAS  II  is 
limited  to  the  treatment  of  single-layer  shells,  however,  a  compromise  shell 
is  considered.  As  outlined  in  Table  3,  the  single-layer  compromise  shell 
possesses  inertial  and  elastic  properties,  and  static-elastic-stability  and 
extensional  elastic-limit  characteristics  that  are  identical  to  those  for  a 
sandwich  shell  with  properties  representative  of  a  stiffened  steel  shell. 

The  two  shells  differ  somewhat  with  respect  to  their  inelastic  flexural  charac¬ 
teristics,  however. 

Figures  2-15  show  response  results  pertaining  to  excitation  of  the  com¬ 
promise  shell  by  a  transverse  plane  wave  of  rectangular  pressure-profile  that 
first  contacts  the  shell  at  9  -  tt  (Figure  1).  The  magnitude  of  the  incident  wave, 

P_,  is  equal  to  four  times  the  shell's  elastic  critical-buckling  pressure,  P  , 
x  c 

which  is,  in  turn,  nearly  equal  to  the  shell’s  axisymmetric  elastic-limit 
pressure,  PQ.  The  width  of  the  rectangular  wave  is  four  shell  radii  (T^  »  4): 
the  incident  wave  loading  is  of  moderately  long  duration.  The  n  •  0-5  circum¬ 
ferential  harmonics  are  considered  and  the  nonlinear  material  behavior  of  the 
shell  is  taken  as  elastic/perfectly  plastic. 

Shown  along  with  inelastic  results  are  results  pertaining  to  nonlinear 
elastic  shell  response,  as  calculated  by  the  present  code  RFP-DYNA  and  the 
small  dynamic  buckling  program  RPF-DBP  [3].  (Good  agreement  between  the 
RPF-DYNA  and  RPF-DBP  results  has  been  a  verification  requirement  for  the 
RPF-DYNA  software.)  Also  shown  for  comparison  purposes  are  linear-elastic 
results.  Finally,  it  is  important  to  remember  that  all  displacement  and 
velocity  histories  are  normalized  to  the  magnitude  of  the  incident  rectangular 

_3 

wave,  P_,  which  is  9.62  x  10  on  a  nondimens ional  basis;  on  a  dimensional 

1  2 
basis,  it  is  this  value  times  pc  . 
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Figure  2  shows  n*0  displacement  histories  with  inelastic  shell  behavior 

both  included  and  excluded.  It  is  seen  that  the  nonlinear-elastic  (N-E)  RPF- 

DYNA  and  RPF-DBP  histories  are  coincident  with  the  linear-elastic  (L-E)  RPF- 

DYNA  history.  When  inelastic  shell  behavior  is  included,  n *  0  response  grows 

rapidly  after  the  n*0  elastic-limit  displacement  w^^-Q.27  P^* -0.0026  is 

reached;  the  motion  slows,  however,  as  the  incident  wave  passes  beyond  the 

shell,  leaving  an  n  ■  0  set  displacement  of  w  *  -6.87  P„  ■  -0.0661.  Note 

os  I 

the  absence  of  oscillatory  behavior  in  all  the  n*0  response  histories;  this 
is  the  result  of  the  heavy  acoustic  damping  provided  by  the  surrounding  fluid 
[1,2]. 


In  [3] ,  the  radial  and  circumferential  displacement  harmonics  w  and  v  , 

n  n 

n  ^  1,  are  expressed  in  terms  of  the  extensional  and  flexural  displacement 

harmonics  e  and  f  as 
n  n 

v  ■  ne  -  —  f 

n  n  n  n  (lf 

w  -  e  +  f 

n  n  n 


Consideration  of  the  en  and  fR  provides  greater  insight  into  the  physical 

processes  at  work  than  does  consideration  of  v  and  w  ;  hence  e  -  and  f  - 

n ..  n  n  n 

histories  are  shown  here. 


Figure  3  shows  n«l  flexural  displacement  histories,  which  constitute 
rigidy-body  motion,  for  both  elastic  and  inelastic  shell  response.  The  coin¬ 
cident  nonlinear-elastic  RPF-DYNA  and  RPF-DBP  histories  follow  closely  the 
linear-elastic  RPF-DYNA  history,  while  the  inelastic  RPF-DYNA  history  deviates 
significantly  from  the  elastic  histories.  This  deviation  is  caused  by 
greatly  increased  n  ■  1  extensional  response  (as  discussed  later) ,  which  feeds 
back  through  the  fluid-structure  interaction  to  affect  rigid-body  motion  [see  the 
last  of  (8)  and  the  last  of  (16)].  The  small  late-time  velocities  exhibited 
by  the  nonlinear  response  histories  are  due  to  the  pressure-gradient  term  in 
the  second  of  (2).  When  this  term  is  omitted,  the  nonlinear-elastic  histories 
coincide  with  the  linear-elastic  history  and  the  nonlinear-inelastic  history 
approaches  an  asymptotic  displacement  of  f^/P^  *  8.0.  Small  effects  of  the 
pressure-gradient  term  are  also  found  in  the  n  £  2  response  histories. 
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Appreciable  differences  among  linear-elastic,  nonlinear-elastic,  and  in¬ 
elastic  histories  are  exhibited  in  Figure  4,  which  shows  n-2  flexural  dis¬ 
placement  response.  A  small  discrepancy  between  the  nonlinear-elastic  RPF- 
DYNA  and  RPF-DBP  histories  also  appears  at  late  times.  For  t  5  7,  inelastic 
f2-response  consists  of  lightly  damped  oscillations  about  a  set  displacement 
whose  value  cannot  be  discerned  from  the  figure. 

Figure  5  shows  n*3  flexural  displacement  histories  for  the  shell.  The 
nonlinear-elastic  RPF-DYNA  and  RPF-DBP  histories  are  in  satisfactory  agreement, 
with  the  linear-elastic  and  inelastic  response  histories  of  comparable  magni¬ 
tude  but  differing  phase.  As  stated  earlier,  the  response  calculations  in¬ 
cluded  the  n*4  and  n  *  5  harmonics;  flexural  displacement  histories  for  these 
are  presented  in  Figures  6  and  7.  The  figures  show  elastic-response  histories 
that,  in  view  of  their  magnitudes,  differ  unimportantly  from  one  another,  and 
inelastic-response  histories  that  oscillate  at  late  times  about  modest  set 
displacements. 

For  elastic  shell  behavior,  it  shown  in  [3]  that  nonaxisymmetric  exten- 
sional  response  is  negligible,  and  the  en  may  be  ignored.  For  inelastic  shell 
motion,  this  is  not  the  case,  as  shown  in  Figure  8.  The  figure  shows  exten- 
sional  displacement  histories  for  n *  1-5  inelastic  shell  response,  as  well  as 
for  n*l  nonlinear-elastic  shell  response;  nonlinear-elastic  shell  response 
for  n  2  would  be  barely  visible  in  this  figiire. 

Figure  9  shows  radial-velocity  histories  at  9  ■  0  and  9  *  it  for  inelastic, 
nonlinear-elastic  and  linear-elastic  response.  It  is  seen  that  the  nonlinear- 
elastic  histories  differ  little  from  their  linear-elastic  counterparts,  but 
that  the  inelastic  histories  deviate  substantially  from  their  elastic  counter¬ 
parts. 

It  is  in  the  strain  histories  of  Figure  10  that  the  substantial  n«l  and 
2  extensional  responses  observed  in  Figure  8  clearly  manifest  themselves. 

They  are  responsible  for  the  largest  set  strains  occurring  on  the  back  side 
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of  the  shell,  with  the  maximum  occurring  at  8  ■  0.  The  proximity  of  the 
inner-fiber  and  outer-fiber  strain  histories  for  each  0-location  illustrate 
the  relative  unimportance  of  flexural  strain  in  the  compromise  shell.  Such 
strain  would  be  even  less  important  in  the  corresponding  sandwich  shell  with 
its  smaller  total  thickness. 


Finally,  Figure  11  shows  deformation  snapshots  at  t  ■  5,  10,  15,  and  20 
for  inelastic,  nonlinear-elastic  and  linear-elastic  shell  response.  These 
are  obtained  by  the  simple  removal  of  rigid-body  shell  motion  during  Fourier 
synthesis,  which  yields  (see  Appendix) 


vd(9,t) 


wd(9,t) 


^[v^t)  +  w^t)}  sin9  + 


n*2 


v  (t)  sin  n9 
n 


wQ(t)  +  lifVj^Ct)  +w1(t)]  cos9  +  wn(t) 

n-2 


cos  n9 


(17) 


The  inelastic-response  snapshots  exhibit  large  axisymmetric  deformation;  the 
modest  nonaxi symmetric  contribution  is  such  that  maximum  deformation  appears 
at  9  »  0  when  t  *  10,  15  and  20.  In  comparison,  the  elastic— response  snap¬ 
shots  exhibit  appreciably  smaller  deformations  dominated  by  ^-displacement. 

3.3  DAA  RESPONSE  CALCULATIONS 


One  of  the  purposes  of  the  present  study  has  been  to  examine  the  accuracy 
of  inelastic  transient-response  solutions  based  upon  approximate  treatment  of 
the  fluid-structure  interaction  in  accordance  with  the  doubly  asymptotic  ap¬ 
proximation  (DAA)  [6,7].  To  this  end.  Figures  12-15  show  RPF-DYNA  and  DAA- 
DYNA  inelastic-response  histories  for  the  shock-wave-excitation  problem  con¬ 
sidered  in  the  previous  subsection. 

Figure  12  shows  n-0  and  n«l  extensional  displacement  histories  for  the 
rectangular-wave-excited  cylindrical  shell.  Large  discrepancies  between  the 
DAA  histories  and  their  RPF  counterparts  are  seen,  with  the  DAA  seriously 
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underestimating  shell  response.  Note  that  the  late-time  asymptote  for  DAA 
n«l  response  is  of  opposite  sign  with  respect  to  its  RPF  counterpart.  The 
impact  of  this  difference  on  strain  response  will  be  seen  shortly.  Large 
discrepancies  between  OAA-based  and  RPF-based  flexural  displacement  histories 
for  n*l  and  2  are  seen  in  Figure  13.  While  these  discrepancies  have  little 
or  no  effect  in  strain  response  comparisons,  they  play  a  substantial  role  in 
kinematic  response  comparisons.  Radial  velocity  histories  at  9-  0  and  9  *>  tt 
are  shown  in  Figure  14.  The  DAA-inelastic  histories  are  seen  to  lie  somewhere 
between  their  RPF-inelastic  and  RPF-elastic  counterparts,  the  latter  appearing 
in  Figure  9.  Hence  the  DAA  calculations  tend  to  underestimate  velocity-response 
magnitudes  on  the  side  of  the  shell  facing  the  incident  wave,  and  to  over¬ 
estimate  them  on  the  back  side  of  t te  shell. 

Finally,  Figure  15  shows  DAA-computed  strain  histories  at  9  »  0,  tt/2  and 
tt.  A  comparison  of  these  histories  with  their  RPF-computed  counterparts  in 
Figure  10  is  rather  startling.  Not  only  is  peak  strain  underestimated  by  a 
factor  of  2.4,  but  the  location  of  peak  strain  occurrence  is  completely 
missed,  the  DAA-calculations  placing  it  at  9  ■  n  and  the  RPF-calculations 
placing  it  at  9  *  0.  This  is  the  direct  result  of  the  sign  difference  in  the 
DAA  and  RPF  late-time  asymptotes  for  n*l  extensional  displacement  response 
(see  Figure  12). 

Figures  12-15  clearly  exhibit  failure  of  the  DAA  as  an  accurate  method 
for  treatment  of  the  fluid-structure  interaction  in  the  present  problem.  The 
reason  for  this  failure  may  be  stated  briefly  as  follows.  The  strain  field 
in  the  shell  is  dominated  by  the  n*0  harmonic.  But  the  DAA,  for  n*0  notion 
of  an  infinite,  circular  cylindrical  body,  is  no  longer  doubly  asymptotic, 
because  the  added  mass  for  such  motion  is  infinite  [17].  Hence  the  DAA  reduces 
to  the  plane  wave  approximation,  which  seriously  attenuates  the  response. 
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Section  4 


CONCLUSION 

This  study  has  involved  the  formulation,  implementation,  and  execution 
of  an  analytical/computational  technique  for  determining  the  geometrically 
and  constitutively  nonlinear  response  of  an  infinite  cylindrical  shell  to  a 
transverse,  transient  acoustic  wave.  The  fluid-structure  interaction  has 
been  treated  in  accordance  with  both  the  exact  residual  potential  formulation 
and  the  doubly  asymptotic  approximation.  Shell  behavior  has  been  treated 
through  utilization  of  the  nonlinear  structural  analyzer  DYNAPLAS  II. 


During  the  study,  the  following  conclusions  have  been  drawn: 

1.  The  formulation  and  implementation  described  in  Section  2  are  accurate 
and  efficient. 

2.  DYNAPLAS  II  is  a  suitable  structural  analyzer  for  problems  of  this  type, 
even  though  it  is  limited  to  the  treatment  of  single-layer  shells; 

DYNAPLAS  III  is  not  suitable. 

3.  For  the  particular  shock-wave-excitation  problem  considered,  axisymmetric 
shell  response  plays  the  dominant  role  regarding  inelastic  shell  behavior. 

4.  In  spite  of  the  dominance  of  axisymmetric  response,  nonaocisymmetric  ex- 
tensional  response  emerges  as  an  important  contributor  to  shell  deforma¬ 
tion;  this  is  in  contrast  to  the  case  of  elastic  shell  behavior,  where 
such  response  is  negligible. 

5.  The  set  profile  produced  in  the  particular  shock-wave— excitation  problem 
considered  resembles  an  oval,  with  the  largest  set  strains  occurring  on 
the  back  side  of  the  shell. 

6.  In  the  particular  shock-wave-excitation  problem  considered,  the  DAA  fails 
as  an  accurate  method  for  treatment  of  the  fluid-structure  interaction. 

7.  Further  studies  are  needed  to  provide  a  broader  assessment  of  the  suit¬ 
ability  of  the  DAA  for  this  class  of  problems. 
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Table  1.  DRY  CHECK  PROBLEM  PARAMETERS 


Shell  of  Figure  1:  PQg  ■  596,  Eq  *  17.5068,  \>Q  ■  0.3,  h  ■  0.05,  a  •  1 

Loading  of  Eq.  (11)  :  p  ml  p  -  1 . 503xl0"4(n2-l) 

o  4  cn 

P  »  10~3  P  -  2.004xl0~7(n2-l) 
n  cn 

tQ  -  8.9,  tl  -  44.5,  tn  -  2Tq 

g  *  386;  P  defined  after  (13);  given  by  (13)  and  (15) 
Calculations  performed  for  0^t<t^+tQ 

The  shell  parameters  shown  pertain  to  a  sceel  sandwich  shell 
with  flexural  stiffness  "El"  twenty-five  times  that  of  a  mono- 
coque  steel  shell  with  h/a  *  0.01  (see  Table  3  for  equivalence 
technique) , 


Table  2.  DRY  CHECK  PROBLEM  RESULTS 

(T  /S  xlO6) 
n  n 


n 

At 

From  (15) 

DYNAPLAS  II 

DYNAPLAS  III 

2 

0.89 

91.9/7.25 

92.6/7.26 

57.9/2.80 

3 

0.89 

32.5/7.69 

33.8/7.69 

4 

0.445 

16.9/7.87 

17.8/7.67 

5 

0.445 

10.5/7.95 

11.7/7.78 
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Table  3.  SINGLE-LAYER  COMPROMISE  AND  SANDWICH  SHELLS 


Single- Layer  Shell 


Sandwich  Shell 


IT 
%  h 

f  1 


_L 

~ h/2 


H  Middle  \*a.y&T  Contributes 

I  No  Mass  or  Stiffness  i 

_L  7TT/  Z2ZZZZZ  l  III  HU  TZZZUZ  Z == 


0.01,  -  -  0.06266 

*  a 


pc  o  pea 

—  -  0.772,  —  »  3.53,  — -  0.02188  7.72,  ~  -  3.53,  -  0.2188 

P  ’  c  ’2  p  c  2 


Inertial  and  Elastic  Properties 


2  2 
p  h  p  c  h  p  h  P  c  h 

’pt"  0.0772,  "Et"  -  0.9620  "pt”  «  —  -  0.0772,  "Et"  -  0.9620 

od  z  pa  t- 

M  pc  a  M  pc  a 


"EI"  «  -  8.017xl0-4 

2  12  a 

pc 


"EI"  -  ^  -  8. 

pc 


017x10 


Static-Elastic-Stability  and  Extensional  Elastic-Limit  Characteristics 


p  c  2  ,  3  - 

P  *  *<“X-T)  (7)  -  2. 405xl0~J 

c  pea 


Pn  c  2  .3  , 

pc  -  25 (—)  (~)  (--)  »  2.405x1 0*J 


.  a  - 

(-)(— ^-)  -  2.462x10"  ** 

a  i. 


N  -  (-)(-4)  -  2. 462xl0~3  ** 
y  a  2 
7  pc 


Inelastic  Flexural  Characteristics 


My  -  4. 103x10 
Ky  -  0.05113 
Mu  -  6.155x10" 


Yield  Moment 
Yield  Curvature 
Ultimate  Moment 


M  -  6.549x10 

y 

Ky  "  0.08170 

M  -  7.099x10“ 
u 


* 


See  Figure  1  for  parameter 
definitions 


19 


von  Mises  yield  condition 
with  v  *  0.3 


NORMALIZED  DISPLACEMENT,  -v*  <t)/P 


a 


Figure  2.  n*0  Displacement  Response  of  Compromise  Shell 

<?I  -  4Pc  *  4Po'  TI  '  4> 
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' l  I  I  I - 1 - 1  I 

0-0  4.0  8.0  12.0  18.0  20.0  24.0  28.0 

NONDIMENS IONAL  TIME,  t 


NORMALIZED  DISPLACEMENT,  f,(t)/P 


Figure  5.  a»3  Flexural  Displacement  Response  of  Compromise  Shell 

«!  *  4fo  *  4V  h  '  4) 
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Figure  6.  n«4  Flexural  Displacement  Response  of  Compromise  Shell 

(PI  -  4Pc  *  4Pq,  Tr  -  4) 
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NONDIMENS IONAL  TIME,  t 


Figure  7.  n-5  Flexural  Displacement  Response  of  Compromise  Shell 

«I  -  4PC  *  4pa'  px  ■  4> 
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NONDIMENS IONAL  TIME,  t 


Figure  S.  Honaxlsynmecrlc  Extanaional  Diap lacement  Response  of  Compro- 
miae  Shell  (?j  •  u?e  *  4Po,  T  -  4) 
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I  . . . 1  ■  ““1 - 1  I  f  I  “ 

0.0  1.0  8.0  12.0  16.0  20.0  21.0  28. 

NONDIMENS IONAL  TIME,  t 


Figure  9. 


RPF-DYNA  Velocity  Response  of  Compromise  Shell 
(Pt  -  «PC  s  4P0,  Tj  ■  4) 
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EXTREME-FIBER  STRAIN 


Figure  10.  RPF-DYNA  Strain  Response  of  Compromise  Shell  (IF  denotes 
inner  fiber;  OF  denotes  outer  fiber;  •  4P£  ••  4Pq,  T 


Deformation  Snapshots  for  the  Compromise  Shell  (actual  deformations  magnified  two  times) 
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Figure  12.  n*Q  and  n«l  Extensional  Displacement  Response  of  Compromise 

Shell  as  Computed  with  RPF-DYNA  and  DAA-DYNA  (PT  ■  4P  S  4P  , 
Tj  -  4)  Ico 
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Figure  13.  n-1  and  n-2  Flexural  Displacement  Response  of  Compromise  Shell  as 

Computed  with  RPF-DYNA  and  DAA-DYNA  (p  -  4Pc  -  4Pq,  ^  -  4) 
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Figure  15.  DAA-DYNA  Strain  Response  of  Compromise  Shell  (IF  denotes 

inner  fiber;  OF  denotes  outer  fiber;  P_  -  4P  ■  4P  , 
Tx  -  4)  I  c  o 


Appendix 

REMOVAL  OF  RIGID-BODY  SHELL  MOTION 


The  deformation  snapshots  of  Figure  11  have  been  produced  through  the 
omission  of  rigid-body  shell  motion  during  Fourier  superposition  of  the  shell 
response  harmonics.  Such  motion  is  given  by  v(9,t)  “-f^Ct)  sin0 ,  w(9,t)  » 
f^(t)  cos9,  where  f^(t)  is  defined  in  (16).  Now  (16)  may  be  inverted  to  yield 


n  ,  .1  . 

■z —  (v  +  —  w  ) 
2..  n  n  n 


n  +1 


(— v  +  nw  ) 
n2+l  Q  n 


(18) 


The  omission  of  f  from  (16)  then  yields 
n 


*  n  .  .1  . 

V  ■  — = —  (v  +  —  w  ) 

n  2 , ,  n  n  n 

n  +1 


*  n  ,  .  1 

wn  "1 —  (vn  +  n  V 


(19) 


n  +1 


which,  for  n  *  1,  lead  to  the  coefficients  of  sin9  and  cos9  in  (17). 
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